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Використання вартащйного принципу найменшо! ди 
до стохастичних систем 


В статье рассмотрена возможность применения принципа наименьшего действия для механических 
систем к стохастическим системам с непрерывным распределением случайной величины. Функция 
Лагранжа выбрана в виде квадратичной формы по функциям распределения Е(х) и плотности 


распределения }(х)=Е(%). В результате получены дифференциальные уравнения, приводящие к четырем, 
хорошо известным в теории вероятности, распределениям: равномерному, линейному, гармоническому и 
экспоненциальному. Результаты получены в случае, когда функция Лагранжа не зависит от случайной 
величины явно, поэтому являются основой для дальнейшего исследования различных форм распределений. 
Ключевые слова: вероятность, функция распределения, плотность распределения, функционал, 
действие, квадратичная форма, стохастические системы, параметр, вариационный принцип. 
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систем з неперервним розподллом випадково! величини. Функця Лагранжа побудована у квадратичнй 
форм! вдносно функщЙ розподлу ЁЕ(х) 1 густини розпомлу ДГ(х)= Е(х). Внасмдок отриман! 
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Лагранжа не залежить в!д випадково! величини неприховано, тому е основою для подальшого досллдження 
разних форм розподрт. 
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Введение 


Существует определенная аналогия в описании динамики механических и стоха- 
стических систем. Динамика механической системы полностью определяется законом 


изменения ее обобщенных координат от времени 4,(#),1 = 1,2,...,8, а начальные 


обобщенные координаты 4;(1,) =4; и обобщенные скорости 4,(1,) = 4» однозначно оп- 
ределяют последующие положения системы. В теории вероятностей такую же функцию 
выполняют функция распределения Е(х) случайной величины х и дифференциальная 
функция распределения }(х)=ЁЕ(х), с той лишь разницей, что начальные условия за- 
меняются краевыми условиями Р(а)=0, Е(Ь) =1 для функции Е(х). Именно эта анало- 
гия может быть положена в установлении уравнений, которым должна удовлетворять 
функция Р(х). 

Для механических систем разработан общий подход их описания, основанный 
на принципе наименьшего действия [1-3] как одном из функциональных методов [4-6], 
используемых в различных областях теоретической физики. 

Целью работы является применение принципа наименьшего действия к сто- 
хастическим системам широкого класса и выяснение общих требований и ограни- 
чений для непрерывно распределенных случайных величин. В частности, установление 
фундаментальных соотношений, которым должны удовлетворять стационарные 
функции распределения и плотности распределения. 

В работе рассмотрены одномерные распределения, но теория легко обобщается 
на многомерный случай независимых случайных событий и позволяет в принципе 
рассматривать зависимые события. 


1 Принцип наименьшего действия в механике 


Известно, что одним из самых общих подходов формулировки закона движения 
механических систем является принцип наименьшего действия [1], согласно которому 
каждая механическая система характеризуется некоторой функцией 1[4(4,4,г) обобщен- 


ных координат 4, обобщенных скоростей 4 = 44/41 и времени 1. Если система имеет 
5 степеней свободы, то а=(4,,4......4,), 9= (4,4,...4.). Функция Ё однозначно оп- 
ределяет закон движения данной механической системы. 

Пусть в моменты времени 1 =[ и #1=5 система имеет координаты а и а"* 


Тогда между а И а система движется так, чтобы интеграл 


5= [`4и49,4,0) (1) 


имел наименьшее возможное значение. Функция Г, называется функцией Лагранжа 
данной системы, а интеграл (1) — действием. 

Рассмотрим систему с одной степенью свободы, так что должна быть определена 
одна функция 4,(1). Пусть 4,(1) искомая функция, для которой 5 имеет экстремум. 


Рассмотрим близкую траекторию 4,(1)+64(1), где 64(1) — вариация функции 4,(1) 
(функция, малая во всем интервале времени от [ до 1, ). Поскольку при Е=Ц и #=1 


все функции 4,(1) +64(1) должны принимать одни и те же значения а и 42) ‚ то 
5, ) = 94.) =0. (2) 
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Изменение 5 при замене 4,(1) на ( 4,(#) +64(1) ) дается первой вариацией 


905 = [ АИХа; +04, 4; +54.1)- |" АИ4,,9,,Г) =| Е: вх РМ , 
Й 1 Е 04. 04, 


Необходимым условием экстремума 5 является условие 05 = 0. Замечая, что 


1 


64 = “54 ‚ проинтегрируем второй член по частям и получим: 
ГА 


В силу условий (2) первый член в этом выражении равен нулю. Остается интеграл, 
который должен быть равен нулю при произвольных ба. Это возможно только в том 
случае, если подынтегральное выражение равно нулю тождественно. Таким образом, 
получаем уравнение 


ЕО. (3) 


При наличии 5 (независимых) степеней свободы в принципе наименьшего дей- 
ствия должны независимо варьироваться 5 функций 4;(1), 1=1.2....,5. Очевидно, что 
получим 5 уравнений вида (3). Эти дифференциальные уравнения в механике на- 
зываются уравнениями Лагранжа. Если функция Лагранжа данной механической 
системы известна, то уравнения (3) {= 1,2,....5 представляют собой уравнения дви- 
жения системы. 

С математической точки зрения уравнения (3) составляют систему 5 урав- 
нений второго порядка для 5 неизвестных функций 4;(1). Общее решение такой 
системы содержит 25 произвольных постоянных. Для их определения необходимо 
знание начальных условий, характеризующих состояние системы в некоторый за- 
данный момент времени, например, задание начальных значений всех координат и 
скоростей а;(#,)=4,„;, 4;(#,) = 4,:, =1.2,...5. 

Рассмотрим некоторые общие свойства функции Лагранжа, которые понадобятся 
в дальнейшем. 

1. Свойство аддитивности функции Лагранжа Г. Пусть механическая система 
состоит из двух частей А и В, каждая из которых замкнута, тогда функции Лагранжа 
частей Сл и Ёв. Если пренебречь взаимодействием между частями, то функция Лагран- 
жа всей системы [= Г, +[,. 

Замечание. Умножение функции Лагранжа системы на произвольный множитель не 
отражается на уравнениях движения. 

2. Функция Лагранжа определена с точностью до прибавления к ней полной 
производной от любой функции координат и времени. Рассмотрим две функции 
Г’(4,4,Р) и Г(4,4,Ё), отличающиеся друг от друга на полную производную по вре- 
мени от какой-либо функции координат и времени /(4,#): 


в у а 
Вычисленные с помощью этих двух функций интегралы (1) связаны соотношением 


й ь 7 А ь о о а 
5'= [‘аила,4,д = [‘аиха, 4+ | 4 Г(а)= 5+ (аль) Гат, 
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т.е. отличаются друг от друга членом, исчезающим при варьировании действия, так 
что условие 95’=0 совпадает с условием 05 =0, и вид уравнений движения 
остается неизменным. 


2 Вариационный принцип в теории вероятностей 


Рассмотрим стохастическую систему с непрерывно распределенной случайной 
АЕ(х) 


величиной х, функцией распределения Р(х) и плотностью {(х)= ‚ которую в 


дальнейшем будем обозначать Е(х). Функция Лагранжа системы Г/Р, Е,х), а действие 
[ " ЦР, Е, 2х. (5) 
Пределы интегрирования могут быть неограниченными х, =-9, х, =+0. 
Произведем варьирование | | =: 5Е+ =: эр =0, учтем 5Е = а, проинтегри- 


руем второй член по частям, получим: 


ов” + | [22 ВР =0. (6) 
ве | 3 (0Е &0дЕ 
В точках х Их, должны выполняться условия 

5Е(х)=6Е(х,)=0. (7) 


В случае неограниченного изменения случайной величины на бесконечности 
функция распределения Р(х) определена значениями ЁР(—5)=0 и Е(+0) =1, поэтому 
Е (—0) =бЕ(+00) =0. 

В силу условий (7) первый член в выражении (6) равен нулю. Остается интеграл, 
который должен быть равен нулю при произвольных значениях 6Р. Это возможно, 
если 

4х ОЕ ОЕ 
При наличии 5 независимых случайных распределений Е, в принципе наимень- 


шего действия должны независимо варьироваться 5 различных функций Ё.(х). Очевидно, 


что получим 5 уравнений вида (8). 
Если функция Лагранжа данной стохастической системы известна, то уравнение 
(8) представляет собой уравнение для определения функции Ё распределения системы. 
В общем случае уравнения (8) для неизвестных Е(х), [=1,2,...5 составляют 


систему 5 уравнений второго порядка. Общее решение такой системы содержит 25 
произвольных постоянных. Для их определения необходимо знание граничных усло- 
вий, например, Е; (—0) =0 и Е(+5) =1, #=1,2,,...5. 


3 Различные виды распределений случаев 
квадратичной формы функции Лагранжа 
Построим функцию Лагранжа по аналогии с «механической» функцией Лагранжа 


2.0 
Г=Т(9)-И(а), где Т = : — обобщенная кинетическая энергия системы, состоящей из 
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одной частицы массой т, И(4) — ее потенциальная энергия во внешнем поле. 
Возможны различные варианты функции Лагранжа. В данной работе ограничимся только 
квадратичной формой переменных К ‚Е. Исследование будем проводить от простого до 
самого общего случая. 

1. Равномерное распределение на отрезке [а,6]. Выберем функцию Лагранжа 
в самом простом виде [=Ё?/2. Как известно, в механике такой выбор функции 
Лагранжа соответствует движению частицы с постоянной скоростью 4=4, =сопя и 
линейному закону движения 4=а1+4, , где 4(0)=а, и 4(0) =а, -— начальные значения 
координаты и скорости частицы в момент времени 1, =0. 

В этом случае ОЁ/ОЁ =0 и уравнение Лагранжа имеет простейший вид ге =0 
Г 

Подставляя ОГ./0Ё = Е ‚ получим Ё =0. Решение уравнения Е =С.х+С, . Используя гра- 


1 
ничные условия ЁР(а)=0, Е(Ь)=1, найдем постоянные интегрирования р. С = г х 
—а —а 
Окончательно, функция распределения Г(х) и плотность Е (х) имеют вид 
х-а : 1 
Е(х) = ‚ 1) =Ео)==—. (9) 
Ь-а Ь-а 


Это хорошо известное в теории вероятностей равномерное распределение или 


распределение с постоянной плотностью =. 
2. Линейное распределение на отрезке [а,65]. Выберем функцию Лагранжа в 
виде [=Ё*/2+2АР , где А — заданный параметр распределения. 
В этом случае ОГ./0Е =2А, 0Г./0Е =ЁЕ и уравнение Лагранжа имеет вид Ё=2А. 
Решение уравнения Р= Ах’ +Сх+С,. Используя граничные условия Р(а)=0, Е(Ь)=1, 
найдем постоянные интегрирования 


1 Аа”) ие) 
— . С, = а. 
Рр-а Рр-а 


С 


Окончательно, функция распределения Г(х) и плотность Е (х) имеют вид 


- , _ р 

Р(х) = Ал? (ав - (ва), (ео) а Ав). 
р-а р-а 

В частности, при А =0 получаем равномерное распределение (9). 


В другом частном случае а=0, 6 =1, А=|! имеем простейшую форму линейного 


(10) 


закона ЁР(х)=х’, /(х) =Ё(х)= 2х. 

3. Гармоническое распределение на отрезке [а,6]. Выберем функцию Лагранжа в 
виде [=Ё°/2-0’Е?/2, где - параметр распределения. Тогда ОГ /0Е = -®°Е, ОГ/ОЁ =Ё, 
получим уравнение осциллятора Е +6’ =0. Решение уравнения ГР = С зшох+ С, созох. 
Используем граничные условия Е (а) =0, Е(Б) =1, найдем постоянные интегрирования 


сз Фа _ соз фа 
О со; орут фа - со; юазт ов зто (а-) 
С зшор + С, с0306=1 п а т а 


С, = р | ео ео. 
со юр зш фа - созфазш ор  зшо(а-Ь) 
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Окончательно, функция распределения Р(х) и плотность Ё(х) имеют вид 


_ зто(а-х) с08%(а-х) 


Е(х) ‚ Го =Ео)=-® 


зпо(а-Ь) эшо(а-Б) 
Замечание. Найденная функция Р(х) нормирована, но еще не является функцией рас- 


‚ (2) >0. (11) 


пределения, поскольку она должна быть неубывающей [7], [8], а это условие никаким 
образом не учтено. Следует потребовать, чтобы изменение параметров ® и а, распре- 
деления (11) удовлетворяли неравенству /(х)>0. 

4. Экспоненциальное распределение на отрезке [а,6]. Выберем функцию Лагранжа 
в виде [=Ё?/2-/2Е?/2, где ^. — параметр распределения. Тогда 0Е./0Е =^.?Е , ОМОЕ=Р.и 
уравнение Лагранжа имеет вид Ё-—^7Е =0. Решение уравнения Е=Се\+Се*. Используя 
граничные условия ЁЕ(а)=0, Е(Ь) =1, найдем постоянные интегрирования 


—> я 
Сие” +С,е № =1 С. = —2 ев "№ — а) 


В 


о +Се м =0 С. = 2е “зв 1 — а) 


Окончательно, функция распределения Р(х) и плотность Ё(х) имеют вид 
т _ 5ПА(х—а) сп (х-а) 
е^ 6-9) 6-9 С (Ь-а)’ 5И^(Ь-а)” 
Для этого распределения имеет место замечание аналогичное, что и для гармони- 
ческого закона (11), именно условие /(х)>0 определяет монотонность функции Е(х). 


Е(х)= = Е(х)=^ 1 >0. (12) 


5. Распределения, основанные на квадратичной форме функции Лагранжа 
(Е,Е) общего вида 
Г=анЕ? +2а» ЕЕ аз +2авЕ + 2аз Ё +азз› (13) 
где а, — элементы симметричной матрицы квадратичной формы есть заданные 
действительные числа. 
В этом случае ОГ/ОЁ =2аЁ+2а,Е +2а., ОЕ10Е =2а,Ё +2а»Е +2а„. и урав- 
нение Лагранжа имеет вид 
а.Ё ана = 0: (14) 
Как видно, часть членов квадратичной формы (13) не входит в уравнение Лагранжа, 


что легко объясняется свойством функции Лагранжа: она определена с точностью до 
прибавления к ней полной производной от любой функции /(Ё,х). 


БВ, Р,з) = Р,Р,з) + “КЕ, 

Запишем измененную функцию Лагранжа, т.е. такую, которая отражается на урав- 

нении Лагранжа , 
Е. =а Е” +а.-Е? +2а3Е. (15) 

Вернемся к уравнению (14). Это уравнение содержит три коэффициента а, а.›,а.з, 
из которых только два независимых и а, *0. Поэтому можно положить а! =1. Тогда 
уравнение (14) примет вид Ё-а„,Е-а.. =0. Из него следуют все ранее рассмотренные 
случаи распределений в зависимости от выбора коэффициентов а.›,а.з [9], [10]. 

Ва, =а., =0 -— равномерное распределение; 

2) а. =0, а, 0 — линейное распределение; 
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3) а, <0, а›, =0 - гармоническое распределение; 

4) а, >0, а, =0 — экспоненциальное распределение. 

Рассмотрим последний случай при а., 0. Составим характеристическое урав- 
нение для уравнения Ё-а„,Е =азз: К? а» =05^? =а,, > М. == а, . Общее решение 
однородного уравнения Ё-а„Е=0 есть Е, =Се”“+С.е“. Частное решение неод- 
нородного уравнения Ё-а»ЁЕ =а.; есть Е»; =—а,./а,›. В результате общее решение 
уравнения 


Ев =Ен+Е5=Се“+Се" аа... (16) 


Пусть случайная величина изменяется на промежутке [0,-+0о), Тогда функция рас- 
пределения должна удовлетворять граничным условиям Е(0)=0, Е (+00) =1. Это оз- 
начает, что в решении (16) следует положить С, =0, С, =-Т, а -а,./а,. =1. В резуль- 
тате, функция распределения Р(х) и плотность ЁР(х) имеют вид 

Ех =1-е^“, Ед =Аем. 
Это хорошо известный экспоненциальный закон распределения. Как видно, он следует из 


принципа наименьшего действия с функцией Лагранжа 
Гы РАКЕ? -2Р). 


4 Пути обобщения теории и граничные условия 


В предыдущем пункте рассмотрен случай квадратичной формы функции Лагранжа 
с постоянными коэффициентами а,. В этом случае теория приводит к дифферен- 


циальным уравнениям второго порядка с постоянными коэффициентами. Такие урав- 
нения хорошо исследованы и их решения не вызывают затруднений [9], [10]. 

Одним из путей обобщения теории является введение переменных коэффициентов 
а,(х). Такое обобщение существенно увеличивает класс функций, относящихся к функ- 
циям распределения Е(х). При этом основное уравнение (14) значительно изменяется, 
например, в уравнении появляется производная Ё(х), которая отсутствовала в прежней 
теории. Кроме того, в уравнении Лагранжа коэффициенты а, становятся частично зави- 
симыми между собой, что представляет значительное ограничение на претендентов на 


функцию распределения. 
6. Распределения, основанные на квадратичной форме функции Лагранжа явно 


зависящей от случайной величины [(ЁР,Ё,х). Рассмотрим функцию Лагранжа в виде 
квадратичной формы, но с переменными коэффициентами 


Г=а (ХЕ? +2а(х)ЕЕ +а›(х)Е? +2а-(х)Е+2а(хЕ+азз(х). 


г 9 .. : 
ЕЕ 241 +212 +2а,Е+2а, -2а»Е —2а.з. 


В результате имеем уравнение 


Г "т р , 
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Это дифференциальное уравнение второго порядка с переменными коэффициентами 
не всегда допускает решение в элементарных функциях и, как правило, требует 
специальных методов решения. 

Для установления единственности решения на отрезке [а,б] обычно достаточно ис- 
пользовать краевые условия Е(а)=0, Е(Ь) =1, которые, в случае неограниченного интер- 


вала изменения случайной величины, заменяются условиями Е(—0) =0, Е(+00) = 1. В неко- 
торых условиях задачи удобно использовать «начальные условия», т.е. Е(—09 =0, ЕСод =0, 
поскольку плотность распределения /(х) = Ё(х) на бесконечности равна нулю. 

Важным свойством функции распределения является ее монотонность. Это 
свойство проще записать через плотность в виде неравенства }(х)>0. Неравенство 
отсекает растущие на бесконечности решения. 

Как видно, существует ряд ограничений на класс функций распределения -— это 
граничные условия (два) для функции распределения и (два) для плотности, а также 
условие монотонности }(х)>0. Эти условия позволяет определять уравнения Лагранжа 
более высокого порядка, чем второй. В этом состоит еще один из путей обобщения 
принципа наименьшего действия. 

Наконец, нами рассмотрена функция Лагранжа в виде квадратичной формы в 
переменных РЁ, Ё ‚ но нет принципиального запрета на более высокие порядки перемен- 


ных РГ, Р. В этом можно рассматривать как еще один путь обобщения рассмотренной 
теории. 


Выводы 


В результате применения принципа наименьшего действия к стохастическим си- 
стемам с непрерывным распределением случайной величины получены уравнения 
Лагранжа для функции распределения случайной величины. Исследована функция 
Лагранжа систем в одномерном случае распределения в рамках квадратичной формы по 
функциям распределения ЁР(х) и плотности распределения /(х)= Е(х) . Функция Лагран- 
жа в приближении квадратичной формы и коэффициентами этой формы, другими сло- 
вами, в отсутствии явной зависимости функции Лагранжа от случайной величины при- 
водит к четырем, хорошо известным в теории вероятности распределениям: равномерному, 
линейному, гармоническому и экспоненциальному. Эти законы распределения по- 
лучены без постановки конкретных физических условий протекания стохастических 
экспериментов. 

Из теоретической механики позаимствована только аналогия: функция распределе- 
ния играет роль «потенциальной энергии» рассматриваемой системы, точнее источника 
внешнего потенциального поля, а плотность распределения, будучи производной от 
функции распределения /(х) = Ё(х) ‚ определяет «кинетическую энергию» системы. 


Функция Лагранжа механической системы определяется как разность между кинети- 
ческой и потенциальной энергиями. Независимость функции Лагранжа стохастической 
системы явно от случайной величины означает, что рассматриваются стохастические си- 
стемы, в которых энергия сохраняется, только переходит из кинетической в потенциаль- 
ную и наоборот. Именно это обстоятельство приводит к четырем фундаментальным 
функциям распределения. Становится очевидным, что остальные известные и неизвест- 
ные законы распределения могут быть получены в рамках данной теории, но только для 
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неконсервативных стохастических систем, т.е. для систем с диссипацией. Это важный 
вывод теории. 
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Г..Р.Мгопепко 
АррПапсе ор те Уапанпопа: Рипсф/е ю Опе-Гтепяопа/ 
эюсйа5йс зуйету 


Ассог4те ю Фе Ше тесВапса| уапайопа[ рипстр/е [1], асйоп Г КЕ, ЕР, х)ах 15 


аррНеЯ №ю зюсвазис зузетз УИ а сопипиои$ 41зе1Бийоп РапсНноп ЁЕ(х) оЁа сазиа| 
уапа Ме х. ТВе Гаотапее Гапсйоп 1$ деегиштей аз а ападгайс Тогт оЁ ЧзлБайоп 
ВипсНоп Е(х) ап @151БаНоп депзиу }(х) = Е (<) 


Г=аЕ” +2а РЕ +а»Е” +2аЁ +2азЕ +азз, (1) 


28 @скусственный интеллект» `2012 


Применение вариационного принципа к одномерным стохастическим системам 6 М 


УПеге а, 1$ зуттейлс таёлх оРа диадгайс Гоги (1). 


Ь о 
Тве Рлпсйопа| | ЦЕ, Е, х)ах 1$ пишитеад Бу Гастапее едиайоп 


в (2) 
4х ОЕ ОЕ 


ш Фе гезой уе оБашед ЧШегепиа| едиайоп$ ог Фе 41561Бийоп Рпсйоп Е, 
\/ысь 1$ заазНед ю Ше Боипаагу соп41опз: Р(а)=0, Е(Ь) =1. 

Тре едчаноп$ Гог Фе а1з61Байоп Рапсйоп Е ‹луе ош уеП-Кпо\уп 4156Бийоп 1а\5 
о? Ше Шеогу оЁ ргоБа Шиу [2-3]. 

1 а, =а. =0 1$ Ше азатбавоп у а соп$апе депзиу 
ия : 1 
о, ==. 
2) а. =0, а. +0 15 Ше Ппеаг Язи йоп 


Е(х) = 


м о 2 
1 о в. 
а р-а 


Е(х) = Ад? + А(аБ - (Б+а)х)+ — 
уПеге А 13а рагатейег оЁ Ше 915 1ийоп. 
3) а, <0, а,. =0 13 Ше Вагмоптс 915 51Бийоп м ИВ а рагатеег © 


_ зшо(а-х) созо (а-х) 


АИ = Е(х)= > 0. 
й зшо(а-р) И) (9) зшо(а-Ь) И) 
4) а,, >0, а,, =0 15 Ше ехропепна[ 915 1БийНоп \ ИВ а рагашеег 2, 
Ре а] сл (х-а) 
ЕР(х) = — = А, 220 
(х) е^ #2) 2.6788) 51). (6 - а) Их) 51). (В — а) Их) 


Рог ехам р[е, ш Ше [аз сазе аЕ а, #0 ап4 хЕ [0,+ю) ме Вауе Ше ме! -Кпо\мп ]а\ 
Ра Е 
Тре гези{$ \уеге оБ{ате4 оп|у Гог сазе жБеп Газгапее псНопз 40 поЕЁ амесНу 


Дереп@ оп гап@ош уама Те. Тве арргоасй 1$ а Баз15 Гог Рагфег Нпдте ап@ сопзисНоп оЁ 
91 егепЕ Роги $ оЁ 4151 фа Ноп$. 
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